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LIFTING  SURFACE  THEORY  FOR  WINGS  IN  LOW  FREQUENCY  SMALL  AMPLITUDE 
YAWING  AND  SIDE  SLIPPING  OSCILLATING  MOTIONS  AT  LOW  SPEEDS 

Huang  Mingke 

(Nanjing  Aeronautics  Institute) 

ABSTRACT 

In  this  paper,  the  mathematical  problem  of 
flow  at  low  speeds  past  a  yawing  oscillating  wing 
is  set  up  by  using  a  coordinate  system  fixed  in  the 
wing.  The  snake-like  tail  vortices  behind  the  wing 
are  investigated  in  detail.  In  the  case  of  low 
frequency  oscillations,  numerical  solution  is  obtained 
by  the  non-steady  vortex-lattice  method  which  can  be 
used  to  calculate  the  rolling  moment  aerodynamic 
derivatives  of  the  wing  due  to  yaw  and  side  slip. 

Some  of  the  computed  results  are  compared  with  the 
experimental  results. 

I.  POSING  OF  THE  MATHEMATICAL  PROBLEM  OF  YAWING 

The  aerodynamic  derivatives  caused  by  yaw  and  side  slip  are 
important  primary  data  for  investigating  the  lateral  stability 
of  an  airplane.  These  derivatives  are  mainly  affected  by  the 
wing.  We  take  a  coordinate  system  oxzy  fixed  with  the  wing  as 
in  Figure  1.  oy  axis  points  upward.  is  the  forward  speed 

of  the  wing,  a  and  are  the  angle  of  attack  and  the  upper  reflex 
angle  respectively.  They  are  both  small  quantities.  We  express 
various  quantities  in  dimensionless  form  with  b,  V^,  b/V^  as  the 
characteristic  length,  characteristic  speed  and  characteristic 
time.  Without  losing  any  generality,  we  shall  only  consider  the 
yaw  oxclllation  of  the  wing  with  angular  velocity  u>  around  the 
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y  axis 


where  t  is  the  dimensionless  time,  w  is  the  circular  frequency. 
We  shall  let  **»■«<  1.  .  Since  the  addition  of  yaw  oscillation  is 
equivalent  to  shifting  the  center  of  rotation,  it  is  not  necess¬ 
ary  to  specifically  study  the  yaw. 


For  low  speed  in  compressible  flow,  corresponding  to  the 
small  perturbation  of  the  moving  coordinate  system,  the  potential 
flow  equation  of  the  absolute  motion  can  be  simplified  as 

yg  .  yg  a**-n 

a*1  ay1  T  a*1  ( i ) 


in  which  4>  is  the  perturbation  velocity  potential.  If  only  the 
first  order  terms  of  a  and  if/  are  retained,  then  the  integral  of 
the  dynamic  equation  of  the  small  perturbation  is  simplified  as 
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(2) 


where  p  and  p  are  respectively  the  pressure  and  density.  The  sub¬ 
script  «  indicates  the  condition  of  free  flow. 

The  boundary  conditions  on  the  xoz  plane  are  as  follows: 

In  the  region  occupied  by  the  wing: 


(3) 


In  the  wake  region,  the  pressure  is  continuous  so  that  Equation 
(2)  must  be  zero.  If  only  the  first  order  term  of  u>  is  kept, 

y  ITlcLX 

then  the  equation  is  simplified  to 


/  a»  _  ag  .  ag 
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(4) 


Here  we  have  already  taken  (dg/a*),„  —  0(vymu)  into  account. 

Outside  the  wing  and  wake  regions: 

Wr-9  -  0  (  5  ) 

At  infinity,  the  perturbation  is  zero. 
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Therefore,  the  problem  of  yaw  oscillation  in  flow  past 
wing  at  low  speed  is  reduced  to  the  solution  of  Equation  (1) 
with  boundary  conditions  (3)  through  (5)  to  be  satisfied.  Com¬ 
pared  with  the  problem  0f  no  yaw,  the  only  difference  is  in  the 
boundary  conditions  of  the  wake  region. 

II.  THE  INSTANTANEOUS  SHAPE  OP  THE  WAKE  VORTEX  LINE  IN  THE 
MOVING  COORDINATE  SYSTEM 


Consider  the  first  order  linear  partial  differential  Equa¬ 
tion  (4).  Take  the  differential (dx,  dz,  dt) in  the(x,  z,  t)space 
ao  that 

dx  —  dt  ^  dx 


—  1  —  1  XU), 


(6) 


Then  Equation  (4)  means  that  d*  »  0  along  this  differential  ele¬ 
ment.  Hence,  it  is  the  linear  element  of  $  =  const  in  the  x,z,t 


space. 


Integrating  Equation  (6)  with  the  initial  conditions  xq, 
we  get 


f*  —  X,  +  0  —  h) 
1*  -  *.  +- 


,,.xatinu><«  otu>tt  _  tfl,m..x  xin  at  _  torm..  cos  t 
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Let  xq,zq  represent  a  point  on  the  rear  edge  of  the  wing,  then 
point  x,z  is  the  position  reached  at  time  t  by  the  vortex 
dragged  out  from  this  edge  point  at  time  tQ.  Eliminating  t 
from  the  equations  in  (7),  we  get  the  instantaneous  equations 
for  the  curve  at  time  t  of  the  snake-like  wake  vortex  line  which 
starts  from  the  rear  edge  point  (xo,zQ).  It  can  also  be  seen 
that  the  instantaneous  shapes  of  the  snake-like  wake  vortex 
lines  from  symmetric  points  on  the  rear  edges  of  the  left  and 
right  wings  (equal  xQ )  are  identical. 


III.  INDUCED  VELOCITY  CAUSED  BY  THE  SNAKE-LIKE  WAKE  VORTEX  LINE 


As  shown  In  Figure  1,  we  shall  calculate  the  induced  velo¬ 


city  vy 


produced  by  a  snake-like  wake  vortex  line  at  a  point 


3 


xl,zl  on  y  *  0  plane  at  time  t.  Assuming  that  the  wake  vor¬ 
tex  circulation  r  does  not  change  with  time,  then  by  using  the 
Biot-Savart  equation,  we  get 


ini 


da  OdL 


Por  symbols  used,  see  Figure  1.  In 
the  above  8  0,  —  A0  +  A0-*  __ 

If  only  the  first  order  terms  are  kept. 


then 
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•vortex  line  direction 
when  oj  =  0 


da  QdL  «•  (  dn0,  —  cos0,A0)dx 
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Figure  1 


Substituting  into  Equation  (8),  we  get 


.  x  r  rf-n-sdn^.x*-** 

->  -  m--  -  r.  11, - h 

r«  ' 


Substituting  the  instantaneous  wake  vortex  line  equation  in  the 
above  and  simplifying,  the  two  integrals  in  Equation  (9)  change 
respectively  to: 


[•  q  -  3dn*,Xx  -«.)  f-iria-Cx-x,) 

r\  ri  v  1  " 

+  -^c««(x  -  x0)  -•—]  +  «.[-  -J-  “  *•> 

- —  dnw(x  —  x»)  +  ——11  dx 

to* 

f  —  cose,  ds  _  (•  _£2l£l  j f -x.cos «(x  -  dn «*  (x  -  x,)l 

},m.,  ri  ri  l  L  " 

+  <i,  [l»da«(x  -  x,)  -  -^j-cos«(x  “  x,)  +  dx 
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where  o,  —  du>,/dt  —  -to,  •aautiatfi 


The  two  integrals  above  apparently  both  converge  but  it  is 
difficult  to  find  their  analytic  forms.  The  interesting  case 
in  airplane  stability  computation  is  low  frequency  oscillation 
but  when  w  -*•  0,  both  integrals  above  are  divergent.  Luckily, 
we  can  separate  out  the  odd  terms  that  make  the  integrals  diverge 
and  discover  that  when  the  induced  velocity  caused  by  the  snake¬ 
like  wake  vortices  (with  equal  circulation  but  opposite  direction) 
at  symmetric  points  on  the  rear  edges  of  the  left  and  right  wings 
is  computed,  these  odd  terms  cancel  out  one  another  and  the 
result  is  still  convergent.  The  derivation  follows  below. 


Take  x  very  large,  eventually  to  approach  infinity,  then 
P 


Consider  first  integral  (10) 


In  the  interval  (xq,  xp  +  xo),  when  w  0,  in  Equation  (10) 

—  £*«»»(*  —  **)  ""  ~T 

*-  — C*  *“  +  0(n?) 

!  v  ;  •  ^ 

-  i 

—  4*-  so*  »(*  —  *0  —  <*»«•(«  ~  **>  +  -^r 

6  -  -M  •  , 


In  the  interval  (x  +  xQ,  »),  since  x  >>  1,  therefore, 

tr  P* 

daft i  **  0,  r|  **  (i  *i? 
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Hence 


This  Integral  can  be  easily  integrated.  When  ai  0,  it  has 
logorithmic  singularity.  Wince  after  such  integration  it 
will  automatically  cancel  out  it  is  not  necessary  to  integrate 
the  explicit  expression.  Equation  (11)  is  simplified  by  the 
same  method.  When  ai  0,  Equation  (9)  simplifies  as: 


5a 


The  integrals  in  the  above  can  all  be  computed.  We  then  let 
Xp  ■+  <*>,  and  include  the  odd  terms  containing  Xp  in  function  f. 
Then  after  lengthy  but  uncomplicated  calculations,  we  finally 
get  the  very  simple  formula 


where 


(13) 


(14) 

(15) 


where 

*;  (16) 

Notice  that  in  Equation  (13)  the  odd  terms  containing  /(x,, 
are  omitted  since  they  will  automatically  cancel  out  when  the 
effect  of  the  free  vortex  with  equal  circulation  but  opposite 
directions  starting  from  symmetric  points  from  the  rear  edges 
of  the  left  and  right  wings  is  computed. 


IV*  VORTEX  LATTICE  METHOD 

Divide  the  wing  into  n  bands  along  the  half  wing  span  and 
then  divide  chordwise  into  m  equal  parts  using  equal  percentage 
line.  The  right  half  wing  is  then  divided  into  N  =  n  x  m  small 
lattices  with  the  left  half  wing  symmetric  to  it.  In  each  lat¬ 
tice  is  placed  a  non-stationary  horseshoe  vortex  as  described 
in  paper  [1],  The  only  difference  is  that  the  section  of  free 
vortices  behind  the  rear  edge  should  be  changed  to  the  snake¬ 
like  vortices  described  in  our  paper,  as  shown  in  Figure  2.  The 
oblique  adjacent  vortex  is  placed  at  the  1/4  chord  point  of  the 
lattice  with  the  control  point  taken  at  the  center  3/4  chord 
point.  The  lattice  span  1  is  a  constant.  According  to  the  com¬ 
putational  formula  of  paper  [1],  it  is  most  convenient  to  take 
b  =  1q/2  as  the  characteristic  length. 

Denote  orf'j  as  a  linear  operator.  It’s  effect  on  the  dimen¬ 
sionless  circulation  r^Ct)  of  the  ith  non-stationary  yaw  oscill¬ 
ation  oblique  horseshoe  vortex  will  yield  the  induced  velocity  in 
the  y  direction  at  its  ith  control  point.  According  to  boundary 
condition  (3),  the  problem  is  to  solve  the  equations 

2N 

T',  —  —a,  —  iign(Ou>„x,<ft  i—  1,2,  •••,2V  (17) 

/- 1 

where  a^x^jZ.^  are  respectively  the  local  angle  of  attack  and 
coordinates  of  the  ith  control  point. 

Expand  the  unknown  circulation  as 

'  r,  -  (r*.),  +  (r>  + +  (7>  +  ( 1 8 ) 

where  the  subscripts  a  and  correspond  respectively  to  the  prob¬ 
lem  of  the  angle  of  attack  and  to  the  problem  of  upper  inversion 
angle.  Apparently ,  (r«..)i  is  the  stationary  dimensionless  cir¬ 
culation  for  the  motion  with  no  yaw. 

Expand  j*'*  as  _  (19) 

urn  i  “  i  t  Jr  , 
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Figure  2 


where  ST t  represents  the  operator  cor¬ 
responding  to  the  non-stationary  oblique 
horseshoe  vortex  with  no  yaw  and  ST* 
is  the  operator  corresponding  to  the 
induced  velocity  increment  caused  by 
the  change  of  the  free  vortices  to  snake¬ 
like  ones  along  x  direction  behind  the 
rear  edge  of  the  wing. 


Obviously, 


Substituting  Equations  (18),  (19)  into  Equation  (17)  and  then 
expanding,  keeping  only  the  first  order  term  of  uymax »  we  then 
get  » 

£  (•h'iCO*-)*)  +  +  .9r\(U>w 

<- 1 

(20) 


■»  —a,  —  i  "  1,  2,  •  •  • ,  2N 


According  to  Section  III, 


(21) 


When  the  matrix  coefficients  l<0«  and  [bjw  are  computed  with 
(14)  to  (16),  we  should  notice  the  fact  that  there  are  two  snake- 
like  free  vortices  with  opposite  directions  in  the  oblique  horse¬ 
shoe  vortex  model. 

•sf/Cr*,)*)  ■■  ~  (**'»■  (22) 

where  («/,)„  may  be  calculated  according  to  Equations  (3.6)  (3.7) 
in  paper  [1]. 


Also  according  to  paper  [1],  when  u>  0 ,  we  have 
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*■•((£»,)  *’)- i  (">*((^)  "■ 

v  l  OwW\ 

In  the  above  equations,  («>,)«  and  are 

with  Equations  (3-7)  and  (3.25)  in  paper  [1]. 


are  calculated 


Substituting  Equations  (21)  to  (24)  in  Equation  (20), 
separating  the  problems  of  the  angle  of  attack  and  that  of  the 
upper  inversion  angle  and  then  separately  equating  the  constant 
terms,  the  sine  terms  (terms  containing  and  cosine 
terms  containing  )°n  the  two  sides,  we  get  the  following 
equations : 

For  problem  of  the  angle  of  attack: 

2  S 

2  “  -4*a« 

J-l 

£  <~MW, -  £  CBMIY.),-  £  [(^f)  2>  -■  2,V( 25 ) 

/-I  J«1  I"1 

The  three  sets  of  linear  algebraic  equations  above  may  be 
solved  one  by  one.  The  first  set  of  equations  correspond  to  the 
problem  of  no  yaw  symmetry.  Thus,  we  only  need  to  solve  for  the 
N  unknowns  on  the  right  wing.  It  may  be  proved  that  the  second 
and  third  set  of  equations  correspond  to  problems  of  anti¬ 
symmetry  and  also  need  only  be  solved  for  the  N  unknowns  on  the 
right  wing. 


Problem  of  upper  inversion  angle: 


2  (»',)<.(*»*  -  -4««gn  O i>*4> 


The  typographic , errors  in  2  arithmetic  signs  of  Equation  (3-25) 
in  paper  [1]  have  been  corrected. 
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These  two  sets  of  equations  are  both  anti-symmetric  and  require 
only  half  wing  solution. 

V.  AERODYNAMIC  LOAD  COMPUTATION 

It  is  not  difficult  to  prove  from  the  pressure  coefficient 
Equation  (2)  that  the  Joukowski  theorem  still  holds  for  non- 
stationary  oblique  horseshoe  vortices.  Consider  the  left  produced 
by  the  jth  vortex  model  of  the  right  wing  as  in  Figure  2. 

The  left  of  BC  vortex  is: 

(Ybc)i  -  +  Vm*,*iXVmbrt)lt 

where  x^z^x^  are  the  coordinates  and  backswept  angle  of  the 

center  point  of  the  vortex.  Substituting  in  the  above  Equation 

(18)  and  keeping  only  the  first  order  term  of  oj  we  get 

ymax 

(  Y«r) i  -  p-i'-C  *  IA>  +  +  tv 

•+•  +  f»(F»M,*j)t F«i(I\»)f ( 27 ) 

where  the  first  term  does  not  contribute  to  the  rolling  moment. 

The  left  of  CD  vortex  Is: 

-  p-FU^Cr,..),  (^^')  +  ^ 

Similarly,  the  left  of  BA  vortex  is 

(Y».)(  -  (A^fL)  +  0(-U.) 

It  may  be  seen  that,  different  from  the  case  of  no  yaw,  with 
yaw,  the  vortices  CD  and  BA  are  also  left-producing  adjacent 
vortices  and  snake-like  free  vortices  which  do  not  produce  lift 
are  only  found  behind  the  rear  edge  of  the  wing.  The  rolling 


moment  around  the  x  axis  may  be  found  by  summing  the  contribu¬ 
tions  of  all  the  vortices. 

VI.  COMPUTATIONAL  RESULTS 


With  the  wing  surface  and  span  as  characteristics,  define 
the  rolling  moment  coefficient  mx*  Again  with  the  half  wing  span 
and  V  as  characteristics,  define  the  dimensionless  time  t  and 
dimensionless  angular  velocity  According  to  the  customary 

method  to  express  aerodynamic  derivatives: 


dm. 


where  8  is  the  slip  angle,  0  —  ip/dt. 


For  many  trapezoidal  backswept  wings,  we  calculate  for  the 
two  center  of  gravity  positions  ir-»0.25  and  0.5  (ratio  of  the 
distance  from  the  front  edge  of  the  average  aerodynamic  chord  bft 
relative  to  the  rotational  center  to  bft).  «!  and  mi  are 
found  by  using  the  transformation  equations  of  the  aerodynamic 
derivatives  when  the  center  of  gravity  is  shifted.  Part  of  the 
numerical  result  and  empirical  estimation  of  the  aerodynamic 
derivatives  are  compared  to  the  experimental  result.  Except 
when  particularly  noticed,  we  use  n  x  m  *  8  x  4  for  calculation 
in  this  paper.  Let  X  be  the  wing  span  chord  ratio,  n  be  the 
root  tip  ratio  and  x  be  the  quarter  chord  line  backswept  angle. 
Figures  3  to  5  are  the  cases  of  plane  wing  with  angle  of  attack 
(ip  *  0),  with  c  as  the  left  coefficient.  Figures  6  and  7  are 

y 

cases  with  upper  inversion  angle  (a  *  0).  Except  for  Figure  7 
where  i|>  Is  measured  in  degrees,  radians  are  used  in  all  other 
cases.  Figure  3  has  been  transformed  to  relative  stationary  axes 
system. 


Relative  body  axes  system  is  used  for  the  rest. 


(e)  X-«\  X  m  2.61,  tj  m  0.23 


<b)  X-«5»,  X  -2.61,  Mr  *  0.25 


Figure  3»  Computational  Figure  4.  Computational 

result  for  *?»/*»  result  for  »■>/«» 

(relatively  stationary  axes) 

- this  paper  - estimation 

paper  [2] 

♦  paper  [3] 

□  paper  [3]  experiment 
<t  paper  [4] 


Key:  (1)  degrees 


Computations  proved  that  for  problems  of  angle  of  attack, 
the  effect  on  the  result  of  whether  or  not  to  consider  the 
change  of  wake  vortices  into  snake-like  is  small,  about  10%. 
However,  consideration  of  this  effect  will  make  the  computational 
result  come  out  closer  to  experiment.  The  time  difference 
derivative  in  problems  of  angle  of  attack  is  entirely  caused  by 
the  change  of  wake  vortices  to  snake-like.  Its  value  is  very 
small . 

For  triangular  wing  with  very  small  span-chord  ratio,  the 
effect  of  the  wake  vortices  on  the  wing  may  be  neglected.  Then 
$  is  equivalent  to  the  theoretical  perturbation  velocity  poten¬ 
tial  of  a  long  thin  body  in  problems  of  stationary  angle  of 
attack  with  no  yaw.  Here  Equation  (2)  becomes  the  theoretical 
pressure  coefficient  equation  of  a  thin  long  body  in  paper  [6]. 
Hence,  our  method  approaches  the  theory  for  a  long  thin  body  for 
this  case. 


Figure  5.  Yaw  derivatives  for  a  plane  wing 

X  *  2.61,  x  *  45°  - this  paper;  o  experiment,  paper 

[4] 
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Figure  6.  Computational  result  for  the  problem  of 
upper  inversion  angle 

X  »  45°,  X  =  2.61,  xT  *  0.25 

- this  paper;  - Weissinger  extended  lift  surface  theory 
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Figure  7.  Comparison  of  computational  and  experimental 
results  for  the  problem  of  upper  inversion  angle 
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VII. 


CONCLUSIONS 


In  this  paper  we  have  only  established  the  lift  surface 
theory  for  wings  in  yaw  oscillation  at  low  speed  and  provided 
the  numerical  method  of  solution  for  the  case  of  low  frequency 
oscillation.  By  separating  the  problem  into  several  symmetric 
and  antisymmetric  problems,  we  have  realized  large  savings  in 
computer  internal  storage  and  in  computer  time. 


The  effect  of  a  side  edge  of  a  wing  changing  into  a  rear 
edge  during  yaw  and  side  slip  has  not  been  considered  in  our 
method,  hence  when  the  effect  of  side  edge  becomes  very  important 
(e.g.,  for  trapezoidal  wing  with  very  small  span-chord  ratio), 
this  method  is  not  applicable. 
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AbdtTMt 

Thia  paper  «cti  up  the  mathematical  problem  o i  Bow  past  a  jawing  oscillating  wing  at 
low  speeds  by  using  a  set  of  axes  fixed  in  the  wing.  The  wake-like  vortices  behind  the  wing 
are  investigated  in  detail.  In  the  case  of  low-frequency  oscillation,  the  problem  is  solved  nu- 
mericaHy  by  the  nonsteady  vortex-lattice  method.  This  method  can  be  applied  to  calculate  rolling- 
moment  derivatives  of  wing  due  to  yaw  and  sideslip.  Some  of  the  results  obtained  are  com¬ 
pared  with  experimental  data. 


SIMPLIFIED  NAVIER-STOKES  EQUATIONS  AND  THEIR  MATHEMATICAL 
PROPERTIES* 
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(Computing  Center,  Academia  Sinica) 

Gao  Zhi 

(Institute  of  Mechanics,  Academia  Sinica) 

ABSTRACT 

In  this  paper,  the  two-dimensional  simplified  Navier- 
Stokes  equations  are  analyzed  qualitatively.  The  authors 
conclude  that  the  above  equations  are  ill-posed  when  the 
tangential  flow  velocity  component  u  <  c(c)  is  the  local 
sound  velocity.  Hence,  instability  arises  in  numerical 
solutions.  In  addition,  the  proper  specification  of 
initial  and  boundary  value  conditions  of  supersonic  vis¬ 
cous  flow  around  a  sphere-cone  body  is  also  discussed. 

Ever  since  1904  when  L.  Prandtl  proposed  the  boundary  layer 
theory,  it  has  greatly  affected  the  research  in  viscous  flow  but 
following  the  development  of  high  speed  machines,  the  boundary  layer 
theory  has  become  insufficient  to  treat  many  important  practical 
problems.  Direct  numerical  solution  of  the  complete  Navier-Stokes 
equations  (abbreviated  as  N-S  equations  below)  has  reaped  fruitful 
result  but  the  volume  of  computation  and  the  amount  of  computation¬ 
al  time  are  still  sizeable.  There  is  a  way  that  has  attracted 
people's  attention,  namely  to  improve  on  the  boundary  equations  or 
to  simplify  the  complete  N-S  equations  in  order  to  reach  the  goal 
of  solving  effectively  practical  problems.  Much  work  has  been  done 

in  this  direction  since  the  60 's  which  may  be  classified  into  three 

aj 

categories:  1.  establishment  of  second  order  boundary  layer  theory; 
2.  thin  shock  wave  layer  approximation,,  and  3.  simplification  of  the 
N-S  equations,  we  would  like  to  point  out  that  when  Re.-*oo  the 

i - 
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thin  shock  wave  layer  equations  cannot  transform  smoothly  into  the 
Euler  equations  while  the  simplified  N-S  equations  will  transform 
smoothly  to  the  non-viscous  equations  for  Re«-*oo  ;  secondly,  the 
simplified  N-S  equations  may  be  solved  with  a  unified  computation¬ 
al  scheme  without  the  need  to  use  the  complicated  method  for  solv¬ 
ing  the  boundary  layer  equations. 

Much  work  has  been  done  in  solving  a  series  of  practical  prob¬ 
lems  with  the  simplified  N-S  equations  and  the  common  problem  en¬ 
countered  is  the  instability  in  the  numerical  solutions  that  has 

c  J3-nj 

appeared  when  the  initial  boundary  value  problem  is  solved.  Many 

A 

people  have  proposed  various  methods  of  avoiding  such  difficulties 
but  few  have  investigated  .k.  their,  utically  the  cause  of  the  instabi¬ 
lity,  the  properties  of  the  ^uations  and  the  proper  specifications 
of  the  boundary  value  conditions .  Our  paper  qualitatively  analyzes 
these  questions  from  the  mathematical  and  mechanical  point  of  view, 
based  on  the  practical  computational  experiences  of  the  authors, 
with  the  purpose  of  discussing  it  with  scholars  both  at  home  and 
abroad. 

1.  SIMPLIFIED  N-S  EQUATIONS 

The  form  of  simplified  N-S  equations  varies  depending  on  how 
the  RHS  is  chosen.  The  practical  difference  is  very  small  when 
the  Reynolds  number  Re^,  is  large.  In  the  following,  we  shall  write 
down  the  simplified  N-S  equations  [27]  derived  by  Gao  Zhi  In  1967 
when  he  analyzed  the  problem  of  viscosity — non-viscous  flow  inter¬ 
ference  (.see  Figure  1). 


-g-  (rp«*)  +  (rttpO  —  0 

P\H  dx  dy  H  dx 

_  \  d  !  d*_\  .  H  du  ,  fi  d  (P\ 

Re.  dy  v  dy  /  Re.  dy  Re.  dy 
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(1.2) 
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*  “  ^/PrRe-,Re.  —  Pr  —  efn/X,a  is  the  blunt  end  stagnation 

point  radius  of  curvature,  k  is  the  object  surface  curvature,  x,y 
are  respectively  the  orthogonal  coordinates  along  and  perpendicu¬ 
lar  to  the  object  surface,  u,v  are  the  corresponding  velocity  com¬ 
ponents.  p->  p->  T%ft%l  represent  respectively  the  pressure,  density, 
temperature,  viscosity  coefficient  and  thermal  conductivity.  Con¬ 
verting  to  dimensionless  quantities,  for  x,y,r  we  use  a;  for  u,v 
we  use  VK  and  for  f .  p,  7\  p.  A  .  we  use  respectively  and 

The  Equations  (1.1)-(1.5)  Include  the  non-viscous  equation 
which  differs  fundamentally  from  the  boundary  layer  equations  in 
its  mathematical  type.  This  will  be  discussed  below. 


2.  TYPE  OP  SIMPLIFIED  N-S  EQUATIONS 


As  is  well  known,  the  classical  boundary  layer  equations  belong 
to  the  parabolic  type  and  are  well-posed  to  be  solved  as  initial 
value  problems.  There  is  no  unified  way  of  specifying  the  simpli¬ 
fied  equations.  Some  people  propose  parabolic  type  equations 
[13-15]  and  some  parabolic-hyperbolic  equations  [16,17].  Paper 
[19]  is  of  the  opinion  that  in  the  vicinity  of  the  object  surface 
(assuming  the  wall  conditions  are  «  ■»  o,  r  -•  0,  T  —  Tw.  )  we  have  the 
property  of  reverse  thermal  conduction.  The  fundamental  type  of  the 
equations  determines  the  way  information  is  propagated  in  the  flow 
field  as  well  as  the  specification  of  the  boundary  conditions  and 
the  numerical  solution,  therefore,  it  is  a  very  important  question. 
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The  simplified  N-S  equations  consist  of  several  first  order 
and  several  second  order  quasilinear  partial  differential  equa¬ 
tions.  It  is  difficult  to  ascertain  its  type  in  the  strict  class¬ 
ical  sense.  Since  the  equations  contain  both  convection  effect 
and  diffusion  effect,  we  can  analyze  the  behavior  of  the  solution 
using  one  of  the  effects  as  the  principle  one. 


In  view  of  the  quasilinearity  of  the  equations,  we  limit  our 
discussion  to  a  small  region  about  a  point.  When  #***  o  ,  we 
introduce  the  following  auxiliary  variable 


r 


l 

Ke. 


C2.1) 


Then  we  can  convert 


(1.1)— (1.5)  to  first  order  equations 


Ox  6y  * 


(2.2) 


where  A,B  are  6x6  matrices,  X,  F  are  6  dimension  column  vectors. 
The  characteristic  equation  of  (2.2)  is 

dec  (ai««  +  af,,')  **  0  *,  j  —  1,  2,  •  •  •,  6  (2.3) 


with  characteristic  root 


(2.4) 


The  result  agrees  with  the  conclusions  of  [16,17]. 


From  (2.4),  it  can  be  seen  that  the  differential  equations 
have  four  multiple  eigenvalues,  reflecting  the  parabolic  property 
with  the  signal  propagating  along  x  =  constant  with  infinite  speed; 
there  are  two  other  non-zero  real  eigenvalues,  representing  the 
hyperbolic  property  with  the  signal  propagating  along  the  character¬ 
istic  directions  with  finite  speed.  Hence,  when  or  when 

diffusion  dominates,  we  may  regard  the  simplified  N-S  equations  as 
a  set  of  equations  with  both  parabolic  and  hyperbolic  properties. 

The  same  conclusion  may  be  reached  with  the  method  of  [25]. 
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On  the  other  hand,  when  y  =  0,  i.e.,  when  (1.1) -(1.5)  are 
reduced  to  first  order  equations,  the  characteristic  roots  are 

it  “  “  Hv/u ,  2m  H  ji  (2.5) 

where  M  is  the  Mach  number.  When  M  >  1,  all  the  characteristic 
roots  are  real  and  the  equations  are  hyperbolic;  when  M  <  1,  X^ 
and  Xjj  are  complex  roots  and  the  equations  are  elliptic.  Hence, 
when  the  diffusion  effect  is  very  small,  Cu  -*■  0),  the  property 
of  the  equations  will  undergo  fundamental  changes.  The  property 
of  the  equations  in  the  intermediate  stages  is  complex  but  we  think 
the  more  important  question  is  whether  the  specification  of  ini¬ 
tial-boundary  values  for  the  Equations  (.1.1 )-(.!.  5)  is  well  posed. 


3.  THE  WELL-POSEDNESS  OP  INITIAL  VALUE  PROBLEMS 


In  [16,17]  it  has  been  pointed  out  the  hyperbolic-parabolic 
property  of  the  simplified  equations  but  from  this  they  reached  the 
conclusion:  it  is  well-posed  to  integrate  downstream  as  an  ini¬ 

tial  boundary  value  problem  from  the  stagnation  point  in  the  flow 
field  of  a  supersonic  flow  around  a  blunt  object.  We  take  a  diff¬ 
erent  view.  We  discover  after  some  analysis  that  when  u  <  c,  the 
Cauchy  problem  is  ill-posed.  For  simplicity,  let  us  consider  the 
major  components  of  the  original  equations  when  the  coefficients 
are  frozen: 


A 


dx  by 


(3.1) 


where  jp  -  (*,  P,  f ,  r) 
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When  u  =  0  and  u  =  c,  the  matrix  A  is  degenerate;  therefore,  it 

is  not  possible  to  explicitly  solve  for  OX  from  Equation  (3.1). 

Qx 

In  order  to  discuss  whether  the  Cauchy  problem  is  well- 

posed  when  x  >  0,  we  notice  the  following  theorem  [20] 


The  necessary  condition  that  the  equations  with  constant 
coefficients 

3*  .  QX  .  n  VX 

Oi  (3.2) 

is  non-degenerate  is  that  the  real  part  of  the  characteristic  values 
of  matrix  B  be  positive.  If  these  characteristic  values  are 
different  from  each  other,  then  the  condition  is  also  sufficient. 


We  use  the  same  method  as  in  [20]  to  derive  that  the  necess¬ 
ary  condition  that  (3.1)  is  non-degenerate  is  that  the  roots  of  its 
characteristic  equation 

<kt(lA  —  C)  —  0 

have  non-negative  real  parts.  The  roots  of  (3-3)  are 


i,  - 1,  -  o 

“  M 
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where  —  (m/c)1  .  It  is  not  difficult  to  see  that  when  u  >  c, 

;  and  when  u  <  c,  i,  <  o,  l,  >  o  •  Therefore,  when  u  <  c, 
(3*1)  do  not  satisfy  the  necessary  condition  for  being  non- 
degenerate.  The  initial  value  problem  is  then  ill-posed. 


The  fact  that  computational  instability  appears  in  the  region 
u  <  c  has  been  mentioned  in  many  works  [13-15 >19].  We  think  that 
the  instability  in  the  difference  solution  is  caused  by  the  ill- 
posedness  of  the  original  partial  differential  equations  when 
u  <  c.  No  instability  appears  in  the  region  u  >  c  for  the  same 
implicit  difference  scheme.  From  a  mathematical  point  of  view, 
the  key  is  that  the  coefficient  matrix  A  of  the  derivative  dX/dx 
in  Equation  (3.1)  becomes  degenerate  in  the  solution  region.  Hence, 
to  overcome  the  instability  in  the  numerical  solution,  we  must 
change  the  structure  of  matrix  A.  For  example,  we  may  treat  Op/Qx 
[13]  explicitly  in  the  x  momentum  equation,  treat  du/dx  [19] 
explicitly  in  the  continuity  equation,  choose  dv/dx  [5]  in  the  y 
momentum  equation  and  use  special  method  [26]  to  prevent  the  devel¬ 
opment  of  perturbation,  etc.  In  fact,  if  we  regard  dp/Ox  in  x 
momentum  equation  or  Qu/Ox  in  the  continuity  equation  as  rela¬ 
tively  fixed,  then  the  characteristic  equation  — C)  —  o 

of  Equation  (3.1)  has  roots 


1> 


«•-*  ■*•-£(! c)>01  ‘•-K**9’*- 


(3.4) 


with  non-negative  real  parts,  hence  satisfying  the  necessary  con¬ 
dition  of  non-degenerate  equations. 


4.  SPECIFICATION  OF  INITIAL-BOUNDARY  CONDITIONS 


As  mentioned  before,  the  simplified  N-S  equations  have  an  ill- 
posed  initial  value  problem  in  the  region  u  <  c.  But  then,  how 
should  the  boundary  conditions  be  specified  so  as  to  guarantee  that 
the  problem  is  well-posed?  .  In  the  following  we  shall  only  quali¬ 
tatively  analyze  the  supersonic  viscous  flow  around  a  blunt  body. 
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In  non-viscous  flow,  the  flow  diagram  is  as  shown  in  Figure  2.  If 
the  initial  value  problem  is  specified  from  the  stagnation  point, 
as  is  done  in  integral  relation  method  [21],  then  we  have  two  prob¬ 
lems.  First,  the  stagnation  point  solution  cannot  be  independently 
established.  It  is  related  to  the  normal  condition  through  the 
singular  line  u  =  c  in  the  flow  field;  second,  it  is  unstable  to 
integrate  downstream  from  the  stagnation  point.  This  instability 
arises  because  the  problem  is  ill-posed. 

For  viscous  flow,  if  the  surface  conditions  are  *  «•  0,  »•»<),  T  *•  Tw 
(no  slip  stream  and  no  injection)  then  the  local  sound  velocity 
line  u  =  c  in  the  flow  field  will  not  intersect  the  object  surface 
but  will  extend  closely  along  the  body  all  the  way  to  the  point  E 
(Figure  3)  near  the  tail  region  at  the  bottom.  In  other  words,  the 
subsonic  region  at  the  top  is  connected  to  the  subsonic  return  flow 
region  at  the  bottom  by  a  very  thin  subsonic  band  near  the  object 
surface.  Will  the  signals  before  and  after  passing  through  this 
band  affect  each  other?  This  is  a  problem  of  much  argument  among 
scholars  In  foreign  countries.  Most  people  think  that  this  effect 
exists  [13-15,19,22],  but  some  argue  from  the  fact  that  according 
to  the  analysis  from  the  point  of  view  of  mechanics  the  signals 
after  will  not  affect  those  before,  hence  the  specification 
initial  boundary  problems  for  the  simplified  N-S  equations  well- 
posed..  ‘'According  to  our  earlier  analysis,  the  former  view  is  more 
reasonable.  As  a  matter  of  fact,  assuming  that  the  simplified  N-S 
equations  are  solved  by  using  the  integral  relation  method,  since 
the  curve  u  =  c  does  not  fall  on  the  object  surface  and  falls  instead 
in  the  bottom  region,  we  may  imagine  that  the  number  N  of  the  strips 
in  the  integral  relation  method  approach  infinity  so  that  the  infin¬ 
ite  many  unknowns  to  be  determined  on  the  stagnation  line  are  asso¬ 
ciated  with  the  Infinite  number  of  regularity  conditions  on  u  =  c. 
When  N  -*■  <*>,  the  quantity  on  the  stagnation  line  x  =  0  is  asso¬ 
ciated  with  the  regularity  condition  on  the  whole  line  u  =  c  which 
means  that  the  "relation"  between  the  signals  before  and  after  does 


exist.  In  principle,  the  solution  must  be  sought  as  a  boundary 
value  problem  in  the  region  A  B  C  D  E  F  A  and  it  is  unreasonable 
to  find  the  solution  along  the  x  direction  by  artificially  fixing 
the  flow  field  at  the  beginning. 

According  to  the  above  analysis,  we  think  the  boundary  condi¬ 
tions  for  solution  in  the  region  A  B  C  D  E  F  A  may  be  given  as 

follows:  Specify  symmetric  condition  x  *  0,  no  condition  on  the 

satisfied 

EF  line,  and  any  additional  conditions  not^at  x«Q  will  be  supple¬ 
mented  by  the  regularity  conditions  of  the  implicit  solutions  on 
curve  u  =  c.  The  y  direction  boundary  conditions  may  be  determined 
by  the  characteristics  of  the  original  partial  differential  equa¬ 
tions.  From  Section  2,  we  know  that  the  equations  have  4  x  =  const 
characteristics,  representing  the  parabolic  property.  Signals  pro¬ 
pagate  along  x  a  const  with  infinite  velocity.  Two  boundary  value 
conditions  must  be  given  on  y  a  0  and  y  =  y  .  Also,  since  the 
characteristic  values  l*(o)  >  o,  A,(o)<0,  therefore,  one  condition  needs 
be  given  at  y  *  0 .  ij(y,)  >  0,  i,(y,)  <  0  ,  therefore,  again  one 

condition  needs  be  given  at  y  =  y  .  If  we  let  the  detached  shock 

s 

wave  of  the  flow  around  the  blunt  body  be  the  unknown,  then  one 
additional  condition  must  be  specified  for  a  total  of  seven  condi¬ 
tions.  When  we  take  three  conditions  on  the  object  surface  and  four 
conditions  on  the  shock  wave,  the  equations  happen  to  be  closed  but 
if  we  keep  the  second  order  derivative  of  v  with  respect  to  y  in 
the  normal  direction  momentum  equation,  then  we  lack  one  boundary 
condition  which  Is  generally  supplemented  by  the  continuity  equa¬ 
tion  on  the  shock  wave  [14]  or  the  normal  momentum  equation  on  the 
wall  surface. 

In  practical  computation,  for  a  narrov;,  long  blut  body,  the 
flow  field  at  the  bottom  and  at  the  top  interact  very  weakly  (.the 
perturbation  gradually  disappears  due  to  viscous  dissipation). 
Therefore,  it  is  possible  to  divide  the  solution  region  into  four 
regions  I,  II,  III,  IV.  I  is  the  blunt  nose  region,  mostly  subsonic 
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II  is  the  body  region  flow  field,  mostly  supersonic  with  the  vicin 
ity  of  the  wall  a  thin  subsonic  region.  Ill  is  the  near-tail 
flow  region  with  a  very  complicated  flow  structure;  IV  is  the  far 
tail  flow  region,  i.e.,  all  supersonic  region.  The  criteria  used 
to  separate  the  four  regions  above  is  that  the  interaction  be¬ 
tween  the  solutions  of  the  two  neighboring  regions  should  be 
negligible.  We  use  the  linear  method  to  find  the  solution  in 
region  I.  Experience  indicates  that  if  the  line  ab  includes  the 
major  subsonic  part,  then  the  effect  of  region  II  on  region  I  may 
be  neglected.  Region  II  then  can  be  solved  as  an  initial  boundary 
problem.  But  there  exists  the  numerical  instability.  A  stable 
numerical  solution  can  be  obtained  simply  by  changing  the  property 
of  the  coefficient  matrix  A  as  pointed  out  in  Section  3>  or  by 
suitably  controlling  the  step  length.  Region  III  should  be  solved 
as  a  boundary  value  problem.  Many  people  still  use  equations  of 
the  type  (1.1)-(1.5)  but  it  is  more  reasonable  to  keep  the  second 
order  derivative  term  along  x  direction.  As  for  region  IV,  it  is 
an  all-supersonic  region.  It  is  reasonable  to  solve  it  as  an  ini¬ 
tial  boundary  value  problem  by  using  Equations  (1.1)  —  (1.5)  -  Cur 
results  for  numerical  solution  of  regions  I  and  II  for  a  sphere- 
cone  composite  object  indicates  that  the  above  method  of  region 
separation  is  workable,  much  more  effective  than  the  method  of 
R.  T.  Davis  [5]  or  time  relation  method.  At  present,  simplified 
N-S  equations  have  been  applied  to  channel  flow  [8],  near  tail 
flow  [9],  far  tail  flow  [10],  compressed  corner  [11],  shock  wave 

attached  surface  layer  perturbation  problem  [12]  and  viscous  flow 

C'3-tSJ 

around  three-dimensional  sharp  and  blunt  nosed  bodies^,  and  the 
effect  is  good. 


Figure  1 
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Abstract 

This  paper  present!  some  qualitative  analysis  tor  simplified  Navier-stnkes  equa¬ 
tions.  Authors  conclude  that  the  system  of  equation*  is  ill-posed  when  tangential  com¬ 
ponent  of  velocity  »  is  smaller  than  e  (e  is  local  sound  velocity)  in  the  flow  field.  There¬ 
fore,  instabilities  have  been  encountered  in  numerical  solutions.  In  addition,  we  also  dis¬ 
cuss  the  problem  of  proper  specification  of  the  initial-boundary  value  for  viscous 
hypersonic  flow  around  a  sphere-cone  body.  •  • 
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EQUATIONS  OP  PLIGHT  MECHANICS  CONSIDERING 
ATMOSPHERIC  TURBULENCE 

Jiang  Yu-xian 


(Beijing  Institute  of  Aeronautics) 


To  study  the  effect  of  wind  on  the  flight  performance, 
the  general  form  of  the  flight  mechanics  equation  with  wind 
must  be  given  [1]  -  [6].  Starting  from  the  difference  between 
the  relative  motion  and  the  absolute  motion,  this  paper  derives 
the  general  form  of  the  flight  mechanics  equations  with  wind. 

Let  G>K,Vt,  W,A,T,M,a>t,at,Q  be  the  vector  forms  of  the 
weight,  air  speed,  ground  speed,  wind  speed,  aerodynamic  force, 
thrust,  aerodynamic  moment,  angular  velocity  relative  to  the 
air,  angular  velocity  relative  to  the  ground  of  the  flying 
vehicle  and  the  angular  velocity  of  the  air  respectively.  The 
body  coordinate  system  is  Oxyz  and  the  inertial  coordinate 
system  is  °gxgygzg'  The  components  of  the  above  vectors  along 
the  coordinates  will  be  indicated  by  the  corresponding  sub¬ 
scripts  <Py  9,  r  are  the  angles  between  the  body  coordinates 
and  the  Inertall  coordinates.  The  cyclic  order  is  <£-*s-*r;a. 
a  is  the  angle  of  attack,  6  is  the  sideslip  angle,  4>  is  the 
drift  angle,  0  Is  the  flight  path  angle,  u  is  the  flow  deviation 
angle,  ct  is  the  sound  speed,  J  is  the  rotary  moment  of  inertia 
(which  components  will  be  indicated  by  the  corresponding  sub¬ 
scripts),  m  is  the  mass  of  the  flying  vehicle,  M  is  the  Mach 
number,  H  is  the  pressure  head,  Q  is  the  flow  rate  of  the  fuel. 
After  certain  simplification,  the  vectorial  form  of  the  flight 
mechanics  equations  with  wind  are 


mdV,/ dt  “  C  4-  A(F„  a,  p,  H ,  M,  •••)  +  T(J,  M,  H,  •••) 
Jdvjdt  —  M(V„  a,  p,  M,  •••) 

v.-v,-w,  a 


(l) 


i 

i 


i 
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The  initial  conditions  are  V-,!,.,  —  V,(o),  «g|,«,  —  «,(0) 

The  dynamic  equation  contains  both  the  absolute  parameters  of 
motion  V  and  w  as  well  as  the  relative  parameters  of  motion 

o  © 

V  and  .  Since 
a  a 

V,-W+V„  w.-fl  +  o,.. 

Hence 

dVJdt  -  (d  VJdt)  +  (dW/it )  +  2  Ox  V”.; 
dujdt  -  (dttjdt )  +  (dO/dt ). 

After  multiplication  with  m  and  J  respectively  and  substituting 
into  (1),  the  following  system  of  equations  are  obtained: 

jy 

m  ~df  “  G  +  A(-V“  “■  w*  M*  •••)  +  T(Q,H,M,  •••) 

-  *  -  mlQ  X  v: 

dt 

dfO  jQ 

^‘5’  “  “«»  M,  H,  •••)  —  / 

V,  —  IV  4*  K,  -  0  4- «. 

The  initial  conditions  are  V.l,^  —  V.(o), —  ».(o) 

In  order  to  solve  the  systems  of  equations  (1)  and  (2),  they 
must  be  transformed  into  their  corresponding  scalar  forms  first. 
To  simplify  the  transformation  of  the  vectors,  the  moving 
coordinate  system  is  taken  to  be  the  same  as  the  body  axis 
coordinate  system  Oxyz.  The  relation  between  V„  and  Oxyz  is 
expressed  by  the  angle  of  attack  « **  —  (0,/(P.)..  ,  and  the 

side-slip  angle  p «» (V4)J(vt),  .  The  relation  between  V 

s 

and  Oxyz  is  expressed  by  the  newly  defined  ascending  angle 

n  •*  (VJr/CV,).,  and  the  lateral  shift  angle  g  (rt),/(vt). 

The  wind  is  a  function  of  time  W(t),  and  is  given  in  the  form 

of  its  three  components  in  the  inertial  system  °gxgygzg>  that 

is,  W  (t),  W  (t),  W  (t).  The  direction  cosines  between 
xg  yg  zg 

the  moving  system  Oxyz  and  the  inertial  system  O-X  y  z  is  given 

o  o  o  o 

in  Table  1. 
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Table  1 


If  the  rotational  motion  of  the  atmosphere  is  not  con¬ 
sidered,  that  is,  JJ=  0,  and  equation  (1)  is  projected  to  Oxyz 
the  dynamic  equations  with  the  absolute  motion  parameters  are 
obtained: 


—  «( 

+ 

G,  + 

+ 

r. 

f"d(Vt)y/dt 

—  I»( 

+ 

G,  +  *4, 

+ 

r. 

md(V,\ldt 

—  m[ 

G,  +  /4, 

+ 

T. 

y»  “  ~ 

dt 

-u.~ 

—  j 

— 

Jcay/dt) 

+ 

M  ,  +  ... 

* 

1 

1 

(/.- 

— 

dw,/dt) 

My  +  .  «  • 

1 

1 

a- 

■ «] 

U  + w. 

+ 

•  •  • 

(3.1) 


The  initial  conditions  are  C^,),( o) •••,«>, —  »,(G)*** 

The  kinematic  equations  are 


(V t)t,  “  C^*)#  cot  Scoter  +  tin  r  «in</r  —  cot  r  tin  ® cot ^>) 

+  (^t)i(«n  r  m9co»<J>  +  cot  r  tin 4>) 

(V $),,  “  (^*)«  tin®  +  (y,),cot  rcoM»  •+•  (P,),(— tinrcot®) 

(?t);  *“  eot®  tin <£)  +  (V,),(  tin  rcot  4  +  cot  r  tin  9  tin  </>) 

+  (vt)i( cot  r  C0t«p  —  tin  r  tin d  tin </• ) 


(3.2) 
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4fe  — ^coiy  — Mttjar  .„,$iay  +  „Jco.r.  -iiferia 9 

co*  r  it  it  it 

(V.)'  -  (F,).  -  (V.),  -  (F#),  -  JF,.  (F.).  -  (F,).  -  »F, 

W ,  “  yyMleo»9eo*<J>  4-  ^**9  —  W/,Ico*dsin</> 

M',  “  lfIf(imr«in^  —  co* y  tin 9  coup)  +  Wftco»rcot9 
4-  JF,((  sin  r  co*<£  4-  co*  y  sin 9  sin </>) 

W,  “  W^C  tin  r  sintfco*^  4-  co*  y  aaip)  —  H',luayco*9 
4-  iF .  4(  co*  r  co*  <p  —  tin  y  sin  3  sin  <p) 
v.  -  4-  (F.)i  4-  (F.)l,  a  *»  -(F.),/(F.)„ 

*~'(Fj./(Fj.,  M-  F./o 

(3.2) 

Similarly,  by  setting  fi  *  0  and  projecting  equations  (2) 
to  the  Oxyz  system,  the  dynamic  equations  with  the  relative 
motion  parameters  are  obtained: 

m  JS-Vih  —  G,  4-  Am  4-  T.  —  i»^  +  «[».(F.),  ~  «y(F.),] 

Jt  dt 

m  —  G,  4-  /f,  4-  T,  —  m  ^jj2-  4-  m[to,(F,),  —  io.(F.),] 

m  ffili  -  G,  4-  A,  +  T,  ~  4-  «[«,(F.).  -  o>,(^.)y] 

dt  dt 

J’d-~  -  -o,  “  A)",".  -  J.,  («."«  ”  4-  •  •  • 

dim,  (  du>  A 

y,  -27  -  -<y.  -  y.)*.«>,  +  y«  -  irj +  M> +  ’  *  * 

y.  ^  -  -(y,  -  yj«olW,  -  y .,(«*,  -  «1)  4-  w.  4-  •  •  • 

dt 

The  initial  conditions  are  (v.),},^  -  (F,,),(o)- «.|«  -  «,(<>)• •  •  • 

The  kinematic  equations  are 

(F,),f  *“  (F.),co*3co*(£  4-  (F,),(  sin  r  »n<£  — •  co*  f  sin 3cos</0 
4*  (F«),(  sinr*»a3co*<£  4-  co»r  *in^») 

(FJ,f  -(U«in9  +  (O,co*rco*»  4-  (F.),(-  *inrco«3) 

(F,)t<  ■“  (F,),(  —  co*3sin^»)  4-  (F,),(  tin  r  co* </>  4*  co*  r  tin  9  tin  <p) 

4-  CF«),(co* rco*<£  —  »nr*in3*in<^) 

J± r  -j£..jia_r  i£_w  rinr  -«,co*r.  421  -  ~  ^  3 

i/  cos  0  dt  it  dt 

v.  -  v'fioi  +  CFJi  +  (F.)i,  «  -  -(f.),/(f.)„ 

P  -  (F.)./(F.)„  M  -  F./« 

(f',),,  -  co.f  +  »\f.  (f,)„  -  w*.  +  <*\>.,  -  c*'.).,  -mf.,  ( *» .  2 ) 
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Let  X,Y,  Z  be  the  drag  force,  the  lift  force  and  the  side 
force,  and  let  $  be  the  rigger’s  angle  of  incidence  of  the 
engine.  Then  the  components  of  A,  T,  G  in  equations  (3.1), 
(4.1)  are  equal.  Prom  the  direction  cosines  relation  between 
X,  Y,  Z,  G,  T  and  Oxyz,  the  following  equations  are  obtained: 


A,  ■»  X  co*  a  cos  ft  +  y  sina  +  Zcosasin0, 

A  ,  •“  X  sin  a  cos  p  +  Y  cos  a  —  Z  sin  a  sin  /? , 

A,  “  —X  sin/?  —  Z cos/?, 

-GsinS,  C,  —  —  G  cos  Seas  r»  G.  «-  G  cos  a  sin  y » 
T,  ~  T  cos  q> , 


Generally,  Ty  *  Tz  =  0.  In  engineering,  the  space  motion  is 
always  simplified  into  the  longitudinal  and  lateral  motions. 

When  studying  the  lateral  motion,  steady  flight  in  the  longi¬ 
tudinal  direction  is  assumed,  that  is,  (V  )  =  (V  )  *  constant, 

St  X  ct  XO 

a  — a,  ,  and  the  elevation  does  not  change.  Simplifying  equa¬ 
tions  (4.1)  and  (4.2)  accordingly,  the  equations  of  motion  in 
the  lateral  direction  with  relative  motion  parameters  are 
obtained: 


/  dio,  dot 


at  at 

(^»)»,  “  C^*)**cosdcos$  4*  (P.).(sinTsin  9 1  cotip  +  cos  y  sin  pi)  ^ 
(^»)«i  ^  cosasin<A)  +  (  Vt) ,  (cosy  cos  <p —  sin  y  sin  9,  sin  </>) 

O',).,  -  O'.).,  +  w,',  (v,)gt  -  (Vg).t  +  Wgt 

dtp  _  ioy  cos  r  dr  dtlt  .  . 

_  —  1  * - ■■  b).  * ■  — -  «n 

cos  a0  dt  it 

V.  -  V +  (y.)>,  p  *,  (y.)./(y.)„,  m  -  ?./« 


(5) 


The  first  equation  in  (5)  can  also  be  expressed  as 
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m^-[G,  +  A.  +  T,—m42*+  mu,  (P*)*]/ (FJh 

Similar  treatment  for  equations  (3.1)  and  (3.2)  gives 
the  equation  of  motion  in  the  lateral  direction  in  terms  of 
the  parameters  of  absolute  motion,  and  the  lateral  shift  angle 
e  is  obtained.  In  general,  fi**t**n  .  Hence  the  drift  angle 
should  be  calculated  according  to  4>  —  <b  —  p  or  <£*»<£  —  e 

Assuming  that  vehicle  is  flying  in  a  fixed  direction, 
that  is,  0  —  ,  and  after  simplifying  equations  (4.1)  and 

(4.2),  the  equations  of  motion  in  the  longitudinal  direction 
expressed  in  terms  of  the  parameters  of  relative  motion  are 


m~~ir 
<X),f  - 

dS/Jt  - 


—  Gs  +  +  T,  —  m  — -  *  -H  mwm(y^)9 

at 

dW,  4™* 

—  G,  +  A,  +  T,  —  m— - /,  ^  - 

{Vt),cotd  —  (K,),« —  (^,Bn9  +  (Vm)y  cos  9  | 
“>*  v.-Vo> Oi  +  (OJ.  «-  U-Vja 


(6) 


The  first  equation  in  (6)  can  also  be  expressed  as 

m  *  "  [G»  +  A* +  T'  “  ■  ^  0\). 

Same  treatment  for  equations  (3*1)  and  (3-2)  gives  the 
equations  of  motion  in  the  lateral  direction  in  terms  of  the 
parameters  of  absolute  motion,  and  the  ascending  angle  n  is 
then  obtained.  In  general,  j?  *  «  ,  and  0  can  only  be  calculated 

according  to  e  — a  —  jj  .  The  initial  conditions  of  equations 
(5)  and  (6)  are  the  same  as  those  of  the  original  equations. 

The  components  of  A,  G,  T  should  be  simplified  correspondingly. 
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The  above  discussion  shows  that  the  relative  motion  and 
the  absolute  motion  of  the  equations  of  flight  mechanics  with 
wind  should  be  distinguished.  The  corresponding  initial 
conditions  should  be  taken. 
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